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ABSTRACT 

We calculate the cosmic microwave background (CMB) bispectrum due to inhomo- 
geneous reionization. We calculate all the terms that can contribute to the bispectrum 
that are products of first order terms on all scales in conformal Newtonian gauge. We 
also correctly account for the de-correlation between the matter density and initial con- 
ditions using perturbation theory up to third order. We find that the bispectrum is of 
local type as expected. For a reasonable model of reionization, in which the universe is 
completely ionized by redshift z r i ~ 8 with optical depth to the last scattering surface 
ro = 0.087 the signal-to-noise ratio for detection of the CMB temperature bispectrum is 
S/N ~ 0.1 and confusion in the estimation of primordial non-Gaussianity is /jvx ~ —0.1. 
For an extreme model with z„ ~ 12.5 and ro = 0.14 we get S/N ~ 0.5 and f^L ~ —0.2. 

Subject headings: cosmic background radiation — cosmology:theory — dark ages, reion- 
ization, first stars — early universe — inflation 
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Introduction 



Secondary anisotropics (jAghanim et al.l 120081 ) in the cosmic microwave background CMB can 
be used to probe the universe after recombination. It is also important to take them into account 
when using CMB to learn about the initial conditions of the universe. One important class of 
secondary anisotropics arises due to the scattering of CMB photons by free electrons during and 
after reionization. In this class, cosmologists have so far concentrated on only one of the terms in 
the second order Boltzmann equation, the product of electron velocity and electron number density 
(v e n e ). I t is known as Sunvaey-Zer dovich effect when the source of electrons is hot gas in galaxy 
clusters (jSunyaev &: Zeldovichlll970l ). If instead of thermal motion veloc ity due to the bulk motion 
of electrons is considered it is known as Ostriker-Vishniac (OV) effect (jOstriker Vishniadll986l ; 
Vishniadll987l ) or kinetic Suny aev-Zel'dovich (kSZ) effect. There are however additional te rms in 
the full second order equations (jBartolo et al.ll2006l . 120071 ; |Pitroull2009l ; ISenatore et al.ll2009al ) which 
also arise due to scattering of CMB photons by electrons and which might be important. Most of 
the work on OV/kSZ effect has focused on the CM B power spectrum . The CMB bispectrum and 
trispectrum were calculated in ICastrol ([2003, 12004 ) , however they calculated the next to leading 
order term which is a six point correlation function of first order terms for the bispectrum. The 
leading order term in bispectrum is a four point correlation function of first order terms. They also 
ignored the de-correlation between the linear and the non-linear quantities in their calculation. 



The CMB bispectrum due to inhomogen eous recombinat i on was calculated in lKhatri &: Wandelt 



(200 9a|) hereafter referred to as KW09 (see ISenatore et al.1 (|2009bl ) for a different approach, also 
Khatri &: Wandeltl (|2009bl )). The same equations need to be solved for inhomogeneous reionization 
and we will follow the treatment in KW09. We will model the inhomogeneous reionization using 
the linear perturbation theory of IZhang et al.1 (J2007J) (hereafter ZHH07). For the recombination 
case the Doppler terms which give rise to OV/kSZ effects were found to be sub-dominant compared 
to the net contribution from the quadrupole and higher order moments of the CMB. We will see 
that this is also the case for reionization. For all calculations the gauge dependent quantities are 
in conformal Newtonian gauge. The cosmological parameters used are baryon density Qb = 0.048, 
cold dark matter density Q c = 0.252, cosmological constant = 0.7, number of massless neutrinos 
N u = 3.04, Hubble constant Hq = 69, present CMB temperature Tqmb = 2.725, primordial Helium 
fraction yn e = 0.24, spectral index of primordial fluctuations n s = 0.95 and as = 0.826. 



2. Inhomogeneous reionization 



We will use the linear perturbation theory of ZHH07 to model reionization. Th e results from 
this model are similar to the bubble model of reionization (jFurlanetto et al.1 120041 ) . Due to the 
fact that reionization is sourced by n on-linear physic s , the validity of any model will have to be 
tested with computer simulations (see lTrac Gnedinl (|2009i ) for a recent review). For our purpose 
the analytical treatment of ZHH07, which captures the essential features of reionization on linear 
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scales, is sufficient. We refer the reader to ZHH07 for details of the model as well as for discussion 
on the validity of this approach. An important input for this theory is a model for the distribution 
of ionizing sources. We will u se the same model used i n ZHH07 which is based on the excursi on set 
treatment of halo formation ([Press Schechterl Il974l ; iBond et al.l Il99ll ; lLacev &: Cold Il993h with 
the minimum mass of a halo given by virial temperature of 10 4 K corresponding to where hydrogen 
line cooling becomes efficient. The spectrum of ionizing radiation is taken to be a power law 



C/3 



(1) 



where [i = \ni/ — lnz^o, v is the photon frequency, v§ = 13.6eV /2irh is the ionization threshold 
for hydrogen, h is the Planck's constant, j(jj>) is the number of ionizing photons emitted at fre- 
quency v per unit parameter [i per collapsed hydrogen atom, /3 is the spectral index of ionizing 
radiation spectrum, £ is the total number of ionizing photons emitted per hydrogen atom and 
Cp = Jq 10 e^ +1 ^dfi is the normalization constant with the spectrum cutoff at /i = 10. We take into 
account Helium reionization by assuming that the first ionization of Helium is identical to that of 
Hydrogen. Although not strictly correct, it should introduce only a small error, unimportant for 
us, since Helium will contribute only about 8% of the total electrons. Second ionization of Helium 



is expected to occur at much lower redshifts (jFurlanetto &: Ohl 120081 ) and will give a negligible 
contribution to the CMB bispectrum. We will consider two different models of reionization arrived 
at by choosing different values of parameter £ in Equation Q] with spectral index /3 = —3. For the 
first model we choose £ = 70 to give the optical depth to the last scattering surface tq = 0.087. For 
the second model we choose £ = 1000 result ing in m = 0. 14 which can be considered a reasonable 



upper limit based on WMAP 5 year results (jKomatsu et al.ll2009f). Figure [T show s the reionization 



history for these two models. We use the RECFAST code (jSeager et al 



1999) to calculate the 



residual mean electron number density after recombination switching to reionization code once the 
electron density due to reionization exceeds the residual value from recombination. The ratio of 
electron number density perturbation to matter density perturbation for comoving wavenumber 
of k = O.OlMpc -1 is plotted in Figure [2j We use the approximate solutions to the perturbation 
equations given in ZHH07 and force the electron bias b e = S e /5 m = 1 once the universe is fully 
reionized. This is a very good approximation to the exact equations of ZHH07 where the bias b e goes 
smoothly to unity. The matter density and hence the electron number density will be non-linear 
on small scales and thus will de-correlate with the linear quantities, for example CMB, on these 
scales. We will use 5 e = b e 5 m on all scales, where b e is calculated using linear theory but S m can be 
non-linear. We will take the de-correlation into account using the third order perturbation theory. 
Note that the non-linearity will be significant only at low redshifts for scales of interest when the 
universe is fully reionized and 6 e = 5 m exactly. Also for the leading term in the bispectrum we 
need to correlate CMB with the electron number density perturbation. However this correlation 
will be small on scales much smaller than the horizon size because CMB traces the perturbations 
at a much higher redshift than that of reionization. Thus the contributions to the bispectrum from 
perturbations in the electron number density will be significant only for scales which are linear and 
where we should expect the linear perturbation theory of reionization to work well. 
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At low redshifts a significant fraction of baryons are expected to be in a diffuse phase called 
warm-hot intergalactic medium (WHIM) (jCen fc Ostrikerlll999l . 12006). The perturbations in these 
baryons are suppressed on small (non-linear) scales compared to the dark matter and for these 
baryons the bias b e should be less than 1. However the contribution to the bispectrum from 
3 < z < 6 is ~ 10% and contribution from z < 3 is ~ few%. This is because of the absence of 
bias b e ~ 10 due to the inhomogeneities in the reionization process once the universe is completely 
reionized and the decreasing optical depth due to the expansion of the universe. At 3 < z < 6 only 
a small percentage of baryons are in WHIM (< 10%). Thus their contribution to the bispectrum 
is less than 1% and the error in assuming b e = 1 negligible. We neglect the ~ few% contribution 
from z < 3 in our numerical calculations. 



3. Cross correlation between matter density and initial conditions 



We will need to calculate the correlation (Px) between the electron number density and linear 
perturbation variables or equivalently the cross-correlation between the non-linear matter density 
(5 m ) and linearly evolved matter density (<5l). We will be interested in the CMB anisotropics 
on scales of angular wavenumber t < 2500, corresponding to the smallest scales predicted to be 
probed by the Planck mission 0. On these scales it is sufficiently accurate to calculate the next 
term in the perturbation expansion which, for Px, means going to the third order in perturbation 
theory for matter density perturbation S m . For Einstein-De Sitter un iverse (Om = = 0) the 



solut i on for 5 m can be written as the following pe rturb ation series (I Vishniac 



1983; iGoroff et al 



I - nt, jj | ita^ — ^ *> ' 

19861 : iMakino et al.lll992l ; Ijain fc Bertschingerlll994l ) (see iBernardeau et al.l (2002j) for a review). 

oo 

MM) = ^Ta^Mk), (2) 



71=1 



where 5\ is the linear matter density perturbation at z = and 5 n is of order 6™, k is the Fourier 
wavenumber and r\ is conformal time. The correlation between the linear and the non-linear matter 
density is then given by 



{a(r])5x(k)6 m (k',r])) = a 2 (77)(<Si(k)<Si(k')> + a^K^k^k')) + higher order terms. 
= (2^) 3 Mk + k') [a 2 (r / )P 11 (fc) + a 4 (7 ? )P 13 (A ; )] , 



(3) 



where -Pn(fc) is the linear matter power spectrum at z = and Pis is the correction given by 



( Suto &; Sasaki 



1991 
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(4) 



: http: //www. rssd. esa. int /Planck 
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Fig. 1. — Reionization history for two models with optical depth to the last scattering surface 
r = 0.087,0.14. 
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Fig. 2. — b e = 5 e /5 m for k = O.OlMpc -1 and optical depth to the last scattering surface tq = 
0.087,0.14. 
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For a general cosmology replacing the scale factor a(n) with the linear growth factor D(rj) in the 



Einstein De-Sitter solution g ives an excellent approximation to the true result (jScoccimarro et al 
19981 ; iBernardeau et al.ll2002l ). Thus we have 



Px(k,v) 



(27r) 3 5 D (k + k')P x (k,rj), 



(5) 



where Pi3(k) is negative signifying de-correlation between the linear and the non-linear density 
fields as expected. 

We plot the ratio of cross power spectrum Px to linear power spectrum Pn n (k, rj) = D 2 (r])Pn (k) 
in Figure [3l For £ < 2500, the CMB bispectrum will get contributions from Fourier modes 
k < 0.4Mpc _1 . It is evident from Figure [3] that for k > O.lMpc -1 the matter density perturbations 
become mildly non-linear (i.e. 0.75 < Px/Plin ^ !)• On these scales comparison with iV-body sim- 
ulations shows that going up to third order in perturbation theory is a very good approximation 
while on smaller scales third order perturbation th eory underestimates the cross-correlation be - 
tween the linear and non-linear matter density fields (jJeong &; Komat su 2006: ICarlson et al.l l2009). 
Taking this de-correlation into account results in replacing the linear power spectrum Pn n in the 
bispectrum expression involving 5 e by the cross power spectrum Px- Equivalently we can define 
an effective transfer function that we can use in the bispectrum expressions derived in KW09. 

Px 



5 e J f = 5, 



lin 



(6) 



4. Bispectrum 

The CMB bispectrum is given by 

B^%lm 3 = (^Lx^^o)^^ (^^o)^^ 2 permutations, (7) 

where a£ m (x,rjo) is the spherical harmonic transform of CMB temperature field and superscript 
indicates the perturbation order. Taking into account all terms that multiply 5 e in the second 
order Boltzmann equation for photons results in the following expression for the angular averaged 
bispectrum (see KW09 for details of the derivation). 



iihi 3 - sr ( £i i2 ^ i R*i4rfs 

m\m2m 3 
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(10) 



where g{rj) is the visibility function, ji is spherical Bessel function, Qjp are the first order CMB 
transfer functions, with @^ = AT^ /T the CMB temperature perturbation, the matrices are 
Wigner 3-jm symbols, 6b = ikVb, \% is the baryon velocity, 6 1 = 3kQ^\ = Q 2 + @po + ©P2> 
@pl are spherical harmonic transform coefficients of the polarization field. P (k) is the power 
spectrum of the initial gravitational potential. 

The Oj, — 6-y term in Equation [10] is the OV/kSZ term which has been the focus of extensive 
research so far. The last term gives negligible contribution. The Yli" ©£"* t erm i n Equation \W\ 
is the new term and it, we will find, dominates over OV/kSZ term. During recombination also 
this term was found to dominate over other terms in KW09. Recently Hernandez-Monteagudo and 
Sunyaev have calculated the effect of this term fo r the scattering of CMB photons in the galaxy 
clusters (|Hernandez-Monteagudo Sunvaevll2010l ). 



5. Numerical Results 



We use CMBFAST (jSeliak Zaldarriagal Il996l ) to calculate all first order quantities. All 



gauge dependent first order quantities are in conformal Newtonian gauge. Figure H] shows B^e 
and contributions from different terms in -Bee • It is clear that X^> 2 ©1 gives the dominant 
contribution. Also the OV term has a sign opposite to that of X^>2 ©£^ t erm - We cutoff the sum 
at t = 1500 which is sufficient for 77 < 7500Mpc. The contribution from 7] > 7500Mpc (z < 3) to 
the bispectrum is small (~few %) because the visibility function is small as evident from Figure 
[5] and also because the perturbation in the ionization fraction is zero since the universe is fully 
ionized by this time and we neglect it. Figure [6] shows the absolute value of bispectrum for our 
two models of reionization for £3 = 200. The bispectrum is clearly of local type. It has however a 
different shape than the primordial bispectrum of local type parameterized by the parameter Jnl- 
The confusion with the estimators of the primordial bispectrum of local type can be quantified by 
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using the following statistic ([Komatsu et al.ll2005l ): 



S y 



E 



prim 



\ prim I 



(11) 



Solving Equation [TT] for /jvl gives the confusion that can be expected if the effects of reionization 
on the bispectrum were ignored. This is plotted in Figure [JJ Planck experiment is expected to have 
error bars on /tvl of ~ 5. The confusion due to the inhomogeneous reionization is much smaller 
and thus can be safely ignored while looking for the primordial non-Gaussianity. We also calculate 
the signal-to-noise (S/N) for the detection of the bispectrum due to inhomogeneous reionization 



(jKomatsu & Spergelll2001 
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(12) 



This is plotted in Figure [8] for our two reionization models and for primordial non-Gaussianity 
with /jvz, = 1. For normal reionization with optical depth tq = 0.087 the S/N for l max = 2500 is 
0.1. For the extreme case with tq = 0.14 we get S/N = 0.5, about 20% more than the primordial 
bispectrum with /jvx = 1. 



6. Conclusions 

We have calculated the leading term in the CMB bispectrum due to inhomogeneous reioniza- 
tion. The bispectrum consists of product of two terms in Equation [HJ B^B^q, i,j = 1,2,3. B$ e 
is due to the correlation of electron number density perturbation with CMB. Bqq is the sum of 
two terms, the correlation of CMB with the peculiar velocity of electrons ( the OV or the kSZ 
term) and the correlation of CMB with all higher order moments of CMB. Since CMB traces the 
perturbations at a much higher redshift, the correlation of CMB with peculiar velocity in Bqq and 
the correlation of CMB with electron number density in B^ & is small on small scales. In particular 
the correlation of CMB with CMB in -B@ e dominates over the peculiar velocity or OV/kSZ term. 
We have found the bispectrum to be of squeezed triangle type, i.e., it peaks where one I mode is 
much smaller than the other two with the contribution to the small I (large scale) mode coming 
from the correlation of S e with CMB and to that of large I (small scale) modes coming from the 
correlation of CMB with CMB. Note that there will be some correlation of CMB with 5 e even on 
small scales due to Thomson scattering. 

If the correlations of CMB with itself are ignored as has been done prior to this work, the 
leading term, which is a four point function of first order terms and which we have calculated, 
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Fig. 3. — Ratio of cross power spectrum between linear and non-linear matter density fields to 
linear matter power spectrum (Px/Plin) for different redshifts. 
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Fig. 4. — B$q and contributions to Bq@ from Yli>2 ®g > OV term 0& — t9 7 titled "Slip" and IlW 
titled "Pol" multiplied by numerical factors as indicated to make them more visible. Clearly the 
dominant contribution to -Bee comes from the Yle>2 ®i term. All quantities are plotted for rj 
near the peak of the visibility function before the universe is fully reionized for both models. Note 
that the sign of OV term is opposite to that of J2e>2 ®£ term. 
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Fig. 5. — Visibility function g{r\) for our models of reionization 




Fig. 6. — Bispectrum for two models of reionization with £3 = 200. The shape of the bispectrum is 
of local type. 
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Fig. 7. — Confusion with the primordial non-Gaussianity of local type parameterized by /tvl (see 
Equation [TTj) as a function of maximum angular wavenumber £rnax • 
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Fig. 8. — S/N for our two reionization models assuming a CMB experiment providing a cosmic 
variance limited measurement of the anisotropics up to £ ma x- Also shown for comparison is the 
S/N for the local type primordial non-Gaussianity with /tvl = 1. 
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would be small. In that case the next to leading order term will be a six-point correlation of only 
the electron number densities and velocities and may be expected to be comparatively important 
since the electron number density an d the y e locity woul d be strongly correlated with each other. 
This six point term was calculated in ICastrd (|2003l . 12004 ) . However in the regime where density is 
slightly non-l i near but velocity is linear, they will get slightly de-correlated. This was ignored in 



Castro! (|2003l . 12004 ) which might have resulted in overestimation of their S/N. They also used an 



instantaneous reionization model which does not include the enhancement in the electron number 
density perturbation due to inhomogeneous reionization, b e , which is expected to be greater than 
one leading to underestimation of their S/N. 

The S/N that we get for the leading term including the correlation s of CMB with itself (Figure 
[8]) is more than an order of magnitude greater than what was found in ICastro ( 2003 . 2004 ) for the 
next to leading order term. It is still below the detection limit of Planck for the models considered 
here. Thus if the reionization occurs at even higher redshifts than our extreme model or if the bias 
b e = S e /5 m is higher than what the model of reionization used by us predicts, then the imprint 
of reionization in the CMB bispectrum may be seen by Planck or post-Planck experiments. We 
would like to point out that there are additional terms in the second order Boltzmann equation, 
the second order electron velocity, CMB monopole and quadrupole, that may also give similar 
magnitude contributions to the reionization bispectrum. The CMB polarization may also get 
important contributions from reionization. However the bispectrum is so small, except in the 
extreme cases, that it is unlikely that these additional terms would change our results significantly. 
More important is the finding that even in extreme cases the confusion with the primordial non- 
Gaussianity of local type is much smaller than one (Figure [7]). Thus inhomogeneous reionization 
should not be a cause of concern when looking for non-Gaussianity in the initial conditions of the 
universe in Planck data. 



We acknowledge correspondence with Jun Zhang on linear perturbation theory of reionization. 
We thank Guilhem Lavaux for checking the cross-correlation between the linear and the non-linear 
matter density perturbations in iV-body simulation. 
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